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We present new long-baseline spectro-interferometric observations of the Herbig Ae star HD 163296 (MWC 275) obtained in the 
H and K bands with the AMBER instrument at the VLTI. The observations cover a range of spatial resolutions between ~3 and 
~12 milliarcseconds, with a spectral resolution of ~30. With a total of 1481 visibilities and 432 closure phases, they represent the 
most comprehensive (u, v) coverage achieved so far for a young star. The circumstellar material is resolved at the sub-AU spatial scale 
and closure phase measurements indicate a small but significant deviation from point- symmetry. We discuss the results assuming that 
the near-infrared excess in HD 163296 is dominated by the emission of a circumstellar disk. A successful fit to the spectral energy 
distribution, near-infrared visibilities and closure phases is found with a model in which a dominant contribution to the H and K band 
emission originates in an optically thin, smooth and point-symmetric region extending from about 0. 1 to 0.45 AU. At a distance of 
0.45 AU from the star, silicates condense, the disk becomes optically thick and develops a puffed-up rim, whose skewed emission can 
account for the non-zero closure phases. We discuss the source of the inner disk emission and tentatively exclude dense molecular gas 
as well as optically thin atomic or ionized gas as its possible origin. We propose instead that the smooth inner emission is produced by 
very refractory grains in a partially cleared region, extending to at least ~0.5 AU. If so, we may be observing the disk of HD 163296 
just before it reaches the transition disk phase. However, we note that the nature of the refractory grains or, in fact, even the possibility 
of any grain surviving at the very high temperatures we require (~ 2100 - 2300 K at 0.1 AU from the star) is unclear and should be 
investigated further. 
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1. Introduction were larger than expecte d by classical accretion disk models 

K> \ dMillan-Gabet et alj|2001l). and were found to be correlated with 

K£ , Herbig AeBe stars (HAeBe) are intermediate-mass young stars, the steUar i uminosity (|Mo^iiE& Millan-Gabej] |20^). This 

j3 . surrounded by large amounts of dust and gas. The distrtbu- supports the idea that the NIR emission is dominated bytbj^ther- 

tion of this circumstellar material remains actively debated. mal emission of hot dust heated by stellar rad i ation . Nat ta"eTaT] 

Various types of models can reproduce the spectral energy dis- {2m) suggested that an inner> optically thin cavity produced by 

tnbution (SED) by considering material in geometrically thin dust sublimation exists inside the disk . At the edge of this region , 

accretion disks ([Hillenbrand et alj [1992J), in a spherical enve- where dufjt condenS ates, the disk is exp ected to puff up be cause 

loEejMiroshmchenko et al. 1997), a puffed-up inner disk rim ofme direct il lumination from the star (Dullemo^SD&OOl; 

(|DullemondetalJ|2001; Isel la & Natta| | 2005|) or a disk plus a [ Isella & Natta 2005)i explaining the size-lumino sity law derived 

halo dVinkovicet al. | | 2006|) . Fitting the SED alone is therefore for Herbig Ae (and late Be) stars by | M onnier & Millan-Gabetl 

highly ambiguous. (|2002l) . Based on a small number of interferometric observa- 

Near-infrai-ed (NIR) long baseline interferometry has al- tionSj simple geometrical models we re proposed to explain the 

lowed us to directly probe the properties of matter within the global morphology of these regions (Qvliilan-Gabet et all l200"lt 

innermost astronomical unit (AU), where key quantities for the lEisner et alJ l2004: Mon nier et al. l l2005L 120061) . However, when 



star-disk-protoplanets interactions are set. The first interferomet- larger set s of data became available, it became clear that the re- 

ric studies of HAeBe showed that the NIR characteristic sizes gions probe d by NIR interferometry are much more complex and 

that a deeper understanding requires the combination of pho- 

Send offprint requests to: benisty@arcetn.astro.it tometric and multi wavelength interferometric measurements at 

* Based on AMBER observations collected at the VLTI (European the mflH-arcsecond resolution and more sophisticated models. 
Southern Observatory, Paranal, Chile) with Arcetri Guaranteed Time 

program 081.C-0124, LAOG Guaranteed Time program 081.C-0794 In this study, we present an analysis of the inner disk sur- 

and open time programs 081.C-0851, 081.C-0098. ~ rounding the HAe star HD 163296 (MWC275). This isolated 



2 



Benisty et al.: The sub-AU disk of the Herbig Ae star HD 163296 



Herbig Ae star is described well by a s pectral type of Al, a 
~30 L luminosity, and a ~2.3 M n mass (Ivan den Ancker et al.l 
19981: iNatta et alj|2004t iMontesinos et al.ll2009l) . It is is located 
at 122^ 1 1 pc and exhibits a NI R excess interpreted as the emis- 
sion from a circumstellar disk dHillenbrand et al.ll992). A l arge- 



scale disk was detected in scattered light (Grady et al 



20001). 



as 



well as at millimeter wavelengths (Mannings & Sargent 1997). 
This inclined disk, traced out to 540 AU, was found to be in 
Kepler ian rotation, and p robably evolving towards a debris disk 
phase dlsella et al.l 120071) . In addition, it also exhibits an asym- 
metric outflow on large scales (>27") perpendicular to the disk, 
with a chain of six H erbig-Haro knots (HH409) that traces th e 
history of mass loss dDevine et all 120001: IWassell et al.l 120061) . 
The spectrophotometric observations that probe the intermediate 
and sm all spatial scales are also compatible with the presence of 
a disk. Douce fet al.l d2006l) studied the warm dust emitting in 
the mid-infrared, located in the surface layers of the interme- 
diate regions of the disk (30-100 AU) and concluded that the 
emission was consistent with a disk that has little flaring. This 
co nclusion is consiste nt with the classification of HD 163296 
bv lMeeus et al.l (12001 1) in their Group II, whose SED can be ex- 
plained by assuming that the inner part of the disk shields the 
outer part from stellar radiation. In the innermost regions, the 
far-UV emission lines have been attributed to optically thin gas 
accreting onto the stellar surfac e, a magnetically confined wind, 
or sho cks at the base of the jet (Deleuil et al. 2005; Swart z etafl 
120051) . Weak X-ray emission (L^/L* ~ 510 6 ) was detected o n 
large scales and attributed to the jet (Giint her & Schm itt 2009). 
In the NIR, the SED time variability was interpreted to be due to 
changes in the inner disk stru cture, on timesca les similar to the 
generation of the HH objects (Sit koet aLl l2008). 



With a disk, signs of accretion, and a bipolar outflow, 
HD 163296 provides an excellent case study to understand how 
circumstellar material is distributed on the sub-AU scale. Its 
NIR disk was resolved by I OTA, PTI, and Keck - Interferometer 
(Mill an-Gabet etal.1 120011: lEisner et all 120091: iMonnier et all 
2005), and at milli -arcsecond (mas) resolut ion with the long 
CHARA baselines dTannirkulam et all 120081) (T08, hereafter). 
T08 found that their observations could not be reproduced using 
models where the majority of the K band emission originates in 
a dust rim, but that an additional NIR emission inside the dust 
sublimation radius could explain the visibilities and the SED. 
They interpreted this additional emission as being produced by 
gas, as sugg ested for other Her big AeB e stars dEisner et alj2 007: 
llsellaet all [20081: iKraus et al.l l2008bl) . The advent of spectro- 
interferometry, as provided by the AMBER instrument at VLTI, 
allows us to simultaneously measure the emission at various NIR 
wavelengths and consequently, to derive temperature profiles for 
the emission, bringing additional constraints on its nature. In this 
paper, we present an observational study of the circumstellar 
disk around HD 163296 at the sub-AU scale, using the largest 
interferometric dataset obtained for a young star so far. The pa- 
per is organized as follows: in Sect. 2, we describe the spectro- 
interferometric observations obtained at AMBER/VLTI and the 
data processing; in Sect. 3, we present the obtained visibilities 
and closure phases. In Sect. 4, we outline a successful disk model 
that reproduces all observables and we discuss its physical ori- 
gin. We summarize our results in Sect. 5. 
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Fig. 1. (u, v) plane coverage of the observations in spatial fre- 
quencies. The observing nights are plotted with different sym- 
bols and the corresponding telescope configurations are reported 
in the figure. 



2. Observations and data reduction 



2.1. Observations at VLTI 

HD 1632 96 was obse rved in the 
strument dPetrov et al 



NIR with the AMBER in- 
the Very Large Telescope 



20071) at .... ... 

Interferometer (VLTI; ' lScholledl20Q7l) . located at Cerro Paranal, 
Chile and operated by the European Southern Observatory 
(ESO). The AMBER instrument allows the simultaneous com- 
bination of three beams in the H and K bands {i.e., from 1.6 
to 2.5 /mi) with spatial filtering. The instrument delivers spec- 
trally dispersed interferometric observables {e.g., visibilities, 
closure phases, differential phases) at spectral resolutions of up 
to 12000. 

In the following, we present observations taken at the low 
spectral resolution mode (R~30) with the 1.8 m Auxiliary 
Telescopes (AT) and the 8.2 m Unit Telescopes (UTs). The data 
were obtained within programs of both guaranteed time and 
open time observations (081.C-0794; 081.C-0098; 081.C-0124; 
081.C-0851). HD 163296 was observed with 14 different base- 
lines of 5 VLTI telescope configurations, during 8 nights from 
May to July 2008. The longest baseline is ~128 m correspond- 
ing to a maximum resolution of 3.5 mas in the K band, and of 
2.7 mas in the H band. In this paper, we use the VLTI nomencla- 
ture to identify the different configurations. A summary of the 
observations can be found in TableQ] including the weather con- 
ditions, the average baseline position angles (Bpa), and projected 
lengths (Bl). The projected baseline is obtained when the vec- 
tor between the two telescopes is projected onto the plane of the 
sky. Because of the Earth rotation, measurements with the same 
physical baseline but at different hour angles correspond to dif- 
ferent projected baselines. 

All the observations were performed with three telescopes, ex- 
cept during the night of the 2008 May 26 when only two 
were available. In addition to HD 163296, three calibrator stars 
(HD 156897, HD 160915, HD 163955) were observed before 
and after each measurement on the scientific target to correct for 
instrumental effects. Their stellar parameters, including their di- 
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Table 1. Log of the interferometric observations. The seeing and tq were measured at 650 nm. 



Date 


Spectral window 


Telescope 


Baseline 


Projected 


Position 


Calibrator 


Seeing 


To 


FINITO 




ljurn] 


configuration 


name 


length B L [m] 


angle B PA [°] 


name 


["] 


[ms] 




24-05-2008 


[1.67-1.80]-[2.06-2.42] 


A0-D0-H0 


AO-DO 


26.7 


81.6 


HD 156897 


1.0 


3.9 


Y 








DO-HO 


53.5 


81.6 
















AO-HO 


80.2 


81.6 










26-05-2008 


[1.62-1.79H2.05-2.31] 


D0-G1-H0 


D0-G1 


69.7 


138.2 


HD 156897 


0.8 


4.4 


Y 


04-06-2008 


[1.60-1.82]-[2.03-2.33] 


E0-G0-H0 


EO-GO 


14.3 


56.5 


HD 156897 


0.6 


3.9 


N 








HO-GO 


28.7 


56.5 
















EO-HO 


43.0 


56.5 










05-06-2008 


[1.61-1.82]-[2.02-2.29] 


E0-G0-H0 


EO-GO 


15.9 


71.1 


HD 156897 


0.8 


4.1 


N 








HO-GO 


31.7 


71.1 


HD 163955 














EO-HO 


47.7 


71.1 










24-06-2008 


[1.65-1.80]-[2.02-2.33] 


U1-U2-U4 


U1-U2 


55.8 


33.7 


HD 156897 


1.1 


1.4 


N 








U2-U4 


82.6 


89.8 
















U1-U4 


122.8 


67.6 










06-07-2008 


[1.63-1.79]-[2.01-2.45] 


D0-G1-H0 


DO-HO 


60.2 


70.4 


HD 160915 


1.0 


2.3 


N 








D0-G1 


69.2 


137.1 
















Gl-HO 


71.4 


7.6 










08-07-2008 


[1.60-1. 82]-[2.01-2.44] 


A0-K0-G1 


A0-G1 


81.4 


121.6 


HD160915 


0.8 


2.9 


N 








Gl-KO 


86.5 


29.9 
















AO-KO 


108.4 


74.9 










19-07-2007 


[1.61-1.80]-[2.06-2.42] 


A0-K0-G1 


A0-G1 


82.4 


124.6 


HD 156897 


1.0 


1.6 


N 








Gl-KO 


88.6 


32.7 
















AO-KO 


115.2 


77.6 











ameters, can be found in Table[2] About 25% of the observations 
were performed using the VLTI fringe-tracker FINITO that uses 
70% of the H band flux to measu re the relative optical pa th dif- 
ference between the light beams (Le Bouq uhTet al.l l2008). 



Table 2. Star and calibrator properties. The latter have been cho- 
sen with SearchCal (http://www.jmmc. fr/searchcal_page. htm) 
and getCal (http://mscweb.ipac.caltech.edu/gcWeb/gcWeb.jsp) 



Star 


V 


K 


H 


Spectral Type 


Diameter [mas] 


HD 163296 


6.9 


4.8 


5.5 


AlVe 


/ 


HD 156897 


4.4 


3.1 


3.1 


F2 


0.8±0.2 


HD 160915 


4.9 


3.8 


3.9 


F5V 


0.7±0.1 


HD 163955 


4.7 


4.5 


4.6 


B9V 


0.3±0.2 



2.2. Photometry 

In addition to this interferometric d ataset, we collected pho- 
tometric data f rom the literature dTannirkulam et all [2008; 
LSitko et al.ll2008h . 

2.3. Data reduction 

The interferometric data reduction was performed following 
Tatul li et all d2007l) . using the amdlib package (release 2.1) 
and the yorick interface provided by the Jean-Marie Mariotti 
Center (JMMC). This led to spectrally dispersed raw visibili- 
ties and closure phases for all exposures of each observing file. 
Not all exposures turned out to provide useful data. In several 
cases, instrumental jitter, insufficient fringe tracking, and unsat- 
isfactory light injection into the instrument led to low contrast 
interferograms of our rather faint source, and we had to select 
the good exposures. Various selection thresholds were examined 
based on the fringe signal-to-noise ratio (SNR) criterion and led 
to the same absolute values for the interferometric observables. 



On the other hand, their accuracy changes with varying selec- 
tions, and the optimal case (i.e., with the smallest errors) was 
obtained with a 20% and 80% best exposure selections for the 
squared visibilities and the closure phases, respectively. In ad- 
dition, data obtained at very high airmass with unstable fringe 
tracking were removed. No selection based on the optical path 
difference (i.e., piston) was performed, since the numbers of use- 
ful exposures could have been a possible source of bias. 
For each night, special care was given to the calibration 
and stability of the AMBER+VLTI instrumental transfer func- 
tion throughout the whole observing period. Measurements of 
HD 163296 were encircled with observations of targets of 
known diameters (see Table [2]) and the transfer function was in- 
terpolated along all calibrations of the night. The errors in the 
calibrated spectral visibilities and closure phases include the sta- 
tistical errors obtained when averaging the individual exposures 
as well as the errors in the calibration star diameters. To account 
for the variation in the transfer function with time, we quadrat- 
ically added the dispersion over the calibrator measurements to 
the errors. The latter is the effect that dominates the error budget. 

With 8 to 14 spectral channels in the K band and 5 to 8 in the 
H band, the total data set (after processing) consists of 1000 and 
304 spectrally dispersed K band visibilities and closure phases, 
respectively, and 481 and 128 H band visibilities and closure 
phases, respectively. The processed data w ill be made ava ilable 
for the community in the OI-FITS format ( Pauls et alj|2005l) on 
the OLBIN websiteQ in January 2010. 



3. Results 

The set of interferometric data presented in this paper is by far 
the largest for a single pre-main sequence star. The correspond- 
ing (u, v) plane coverage is shown in Fig. [TJ In this section, we 
present a summary of the results and describe the main charac- 
teristics of the inner region of HD 163296 as measured by NIR 



1 http://olbin.jpl.nasa.gov 
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Fig. 2. Squared visibilities against spatial frequencies B/A in the K band (left panel) and the H band (right panel). Different symbols 
correspond to different configurations, that are reported in the upper right corner of the left panel. 
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Fig. 3. Broad-band squared visibilities against spatial frequencies B/A in the K band (left panel) and the H band (right panel). 
Different symbols and colors correspond to different configurations that are given in Fig. [2] The IOTA points are overplotted in the 
right panel (empty circles). 



interferometry. For the sake of clarity, we separate the H and K 
band results in most of the figures below. 

3.1. Visibilities 

Figure [2] presents the spectral visibilities as a function of the 
spatial frequency (i.e., the ratio of the projected baseline length 
to the wavelength of the observation). Figure [3] similarly shows 
broad-band visibilities, obtained when averaging over all spec- 
tral channels, for the K band (left) and the H band (right). One 



can immediately see from the figures that at these spatial reso- 
lutions, the circumstellar matter around HD 163296 is resolved. 
Within the error bars, the H band visibilities vary with baseline 
qualitatively in the same wau as those of the K band. However, 
at spatial frequencies higher than 20 m///m, the H band V 2 are at 
least 25% higher than those in the K band. 

The errors in individual points vary significantly from night 
to night. However, some of the scatter is not caused by the un- 
certainties in individual measurements, since different observa- 
tions were performed with baselines of similar projected base- 
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lines (equivalently, at similar spatial frequencies) but different 
position angles, therefore sampling distinct orientations on the 
sky. This effect means that the source geometry deviates from 
a circular one, and can instead present an elongated shape, as 
expected for an inclined disk. 

We note that in the H band, at the shortest spatial frequencies 
(below 10 m///m), the AMBER visibilities do not reach unity. 
This can be due to fl ux from an extended halo as suggested by 
iMonnier et al.l d2006l) from their IOTA H band measurements. 
All the IOTA points, except one, are consistent with our mea- 
surements (see Fig. [3] right). On the other hand, this effect is not 
clearly seen in the K band, where visibilities are close to unity 
at short baselines {e.g., V 2 =0.94+0.03 at a 13.8 m baseline). 

The observations show that V 2 also depends on wavelength 
(Fig. [2]). This can be related to the physical extension of the emit- 
ting region at various wavelengths and the existence of tempera- 
ture gradients within it. Since the interferometer resolution also 
changes with wavelength, it is natural to visualize this depen- 
dence using a geometrical model to convert the measurements 
in angular sizes while taking this effect into account. With this 
aim, we fit the V 2 in each spectral channel for each measure- 
ment, using a ring of uniform brightness (with a 20% thickness). 
Although the circumstellar material is mainly responsible for 
the emission in the H and K bands, the star also contributes to 
the measured fluxes and visibilities. We estimate a contribution 
from the stellar photosphere of ~33% in H and -14% in K (see 
Sect. |4.1| for details). The flux from the circumstellar matter was 
then computed for each AMBER channel as the difference be- 
tween the observed NIR flux (T08) and the photospheric flux. 
Figure |4] gives an example of the wavelength dependence of the 
size obtained on the three baselines of the A0-D0-H0 configu- 
ration, in the two extremes cases (observed maximum and mini- 
mum variations with wavelength). 

For the majority of the measurements (except those obtained 
at very small baselines with the E0-G0-H0 configuration), the 
size of the circumstellar matter slightly increases with wave- 
length. To quantify the wavelength dependence of the visibil- 
ity, we studied its variation over only the K band, H band, and 
over both bands together. Across the K and H band separately, 
24% and 31%, respectively, of the measurements show a chro- 
maticity above the 2-cr level. Over the entire spectral range, a 
stronger effect is expected due to the greater wavelength inter- 
val, 62% of the data show a chromaticity greater than the 2-cr 
level. This is consistent with the case shown in Fig. |4] where the 
characteristic size of the emitting region increases between the 
H and K bands. This trend is also observed in the spectrally dis- 
persed measurements acr oss the Br T emissio n line obtained with 
the Keck-Interferometer (Eisner et al . 200§. 

3.2. Closure phases 

The closure phase is a quantity that can be derived from interfer- 
ometric observations with at least 3 telescopes. By combining up 
the phases of the fringes obtained with 3 telescopes, the atmo- 
spheric disturbances are cancelled out. Consequently, the sum 
of the three phases (the closure phase) is atmosphere-free, i.e., 
independent of the phase fluctuations. It is related to the degree 
of asymmetry of the observed brightness distribution: a point- 
symmetric object will have zero closure phases, while a non-zero 
measurement is indicative of a deviation from point-symmetry. 
The sign of the CP is derived from the way that it is calculated 
(clockwise or counterclockwise). Physically, the CP depends on 
the directions that are sampled by the individual baselines of 
the configuration since they probe different asymmetries of the 
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Fig. 4. Example of wavelength dependence for the squared visi- 
bilities obtained on the three baselines during observations with 
the A0-D0-H0 configuration. In each panel, the characteristic 
size of the emission is given at each wavelength across the H and 
K band. These sizes are derived from a ring model of uniform 
brightness to account for the change of resolution with wave- 
length. The triangles show the maximum chromaticity over the 
whole wavelength range, while the circles give the minimum 
variation. For each measurement, the corresponding projected 
baseline length and position angle are indicated in the top right 
corner. 



brightness distribution, which in turn strongly depend on the ra- 
tio of the star to circumstellar matter flux. 

One or more closure phase (CP) measurements were ob- 
tained for each observing night, except for the night of 2008 
May 26. In total, we obtained 304 CP in the K band, and 128 
in the H band. This is an enormous improvement over the ex- 
isting datasets for a young star. The full set of measured CP 
across the K and H bands is shown in Fig. [5] plotted against 
the maximum spatial frequency B max //1, where B max refers to 
the projected length of the longest baseline in the corresponding 
configuration. 

The CP signal depends on the spatial resolution achieved by 
the interferometer (unresolved sources are centro-symmetric). 
Figure [6] show the broad-band CP in K and H plotted against 
the maximum spatial resolution A/B lmx , achieved with the cor- 
responding three-telescope configuration. It can be seen that the 
level of the CP signal increases with the power of resolution 
(see, for example, the configurations of aligned telescopes along 
the same B PA , E0-G0-H0, and A0-D0-H0). Table g] gives the 
corresponding broad-band CP averaged over all the measure- 
ments as well as the achieved maximum spatial resolution for 
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Fig. 5. Measured closure phases for each configuration (reported in the top left corners) plotted versus maximum spatial frequency. 
In all panels, except for U1-U2-U4 for which only K band data were useful, the measurements include both the H and K bands. 



each telescope configuration. Since we checked that all the CP 
were calculated in the same way, the relative signs between mea- 
surements obtained with various configurations are relevant and 
should be reproduced by model fitting. 

The CP is close to zero for the short linear array (E0-G0-H0): 
we do not detect any meaningful deviations from zero {i.e., any 
flux asymmetries) at these low spatial res olutions (10.8 and 8.5 
mas in the K and H band, respectively). iMonnier et alj (2006) 
detected CP signals below < 5° for 12 Herbig AeBe stars with 
the IOTA-3T interferometer, including HD 163296 for which 
they derived 0.6+0.4° at a resolution of 11.8 mas. Our obser- 
vations at similar resolution provide comparable results. 

At higher resolutions, we found that the CP signal is not zero, 
i.e., that the emission is no longer centro-symmetric. However, 
the departure of the closure phase signals from zero is small, 
with a maximum average CP in K band of 1 1.7+2.9° and in H 
band of -5.8±2.5° (see Table[3]). As we discuss in the following 
section, this level of asymmetry is not compatible with strongly 
skewed distributions for the circumstellar material in the inner- 
most regions surrounding the star. 

Except in one case (A0-K0-G1), the CP does not vary much 
with varying hour angle (see Fig.fTZl). This means that, consid- 
ering the change in maximum resolution and position angle that 
occurs when varying the hour angle {i.e., when changing the di- 
rection and the projected baseline length Bl), the level of asym- 
metry does not change much. Consequently, the circumstellar 
matter must have a rather smooth azimuthal brightness distribu- 
tion. 



The H band CP are slightly lower than those measured in the 
K band, as expected since the emission in H probably originates 
in a more compact region, as most of our data suggest, than the 
one emitting the K band flux. However, considering the large er- 
ror bars, this effect is hardly significant. Over the K band and H 
band separately, 12% and 15% of the CP, respectively, show a 
variation above the 1-cr level and all measurements are consis- 
tent with variations within the 2-cr level. Considering the entire 
range of wavelength {H and K together), 45% of the CP show a 
variation with wavelength above the 2-cr level. 

4. Modelling the interferometry results 

The analysis of visibilities and closure phases requires the as- 
sumption of a model for the brightness distribution on the plane 
of the sky to be compared with the observations. Even when the 
comparison is successful, it is impossible to know if the solu- 
tion is unique, and because models are difficult to compute, it is 
practically impossible to explore all possibilities. Here, we are 
guided by the current paradigm that the NIR emission in excess 
of the photospheric one is produced by the inner parts of a cir- 
cumstellar disk. 

4.1. Stellar parameters 

To model the interferometry results, one needs to know the unre- 
solved contribution of the star to the total flux at the wavelength 
of interest. We compute it from the observed magnitude in the 
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Fig. 6. Broad-band closure phases plotted against the maximum spatial resolution achieved for all telescope configurations (K band, 
left; H band, right). The configurations and the corresponding symbols are reported at the bottom of the left panel. 



Table 3. Values of the closure phase as averaged over all the 
measurements obtained with the various configurations. The 
dates of the observations, the wavelength range, and the achieved 
maximum spatial resolution are reported. 



Configuration 


Date 


Band 


Resolution 
[mas] 


CP 

o 


A0D0H0 


24/05/08 


K 


7.9 


6.5±1.4 






H 


5.8 


4.3±1.8 


E0G0H0 


04/06/08 


K 


11.6 


0.0±2.1 






H 


9.0 


-1.5+1.3 




05/06/08 


K 


10.0 


-2.0±1.7 






H 


7.9 


0.2+1.1 


U1U2U4 


24/06/08 


K 


3.8 


-0.5+10.1 


D0H0G1 


06/07/08 


K 


7.0 


-3.3+1.6 






H 


5.1 


-5.8±2.5 


A0K0G1 


08/07/08 


K 


3.3 


11.7+2.9 






H 


3.2 


0.3+4.6 




19/07/08 


K 


4.1 


3.2+2.6 






H 


3.0 


3.3±3.7 



I band (i.e., 6.71 + 0.02 6, with very little variability over a pe- 
riod o f about 20 years (Ide Winter et alJBOOU iTannirkulam et alj 
2008), assuming ZAMS colors for a Al star and an extinction 
Ay = 0.25. The resulting stellar fluxes are about 1.4 Jy in K and 
2.2 Jy in H. The observ ed NIR fluxes display a moderate vari- 
ability dSitko et al.l l2008), so that the stellar contribution may be 
in the interval 14-18% in K and 33-37% in H, respectively. The 
effect of variability on the interferometry results and the desir- 
abilit y of performin g simultaneous photometry were discussed 
bv lSitko et al.l (120081) . This, unfortunately, is practically impos- 
sible with AMBER/VLTI, and we do not know the values of the 
total flux at the time of the observations. In this paper, we adopt a 
stellar contribution of ~14% and ~33% to the observed K and H 
band fluxes, respectively, i.e., on the lower side of the estimates. 
These numbers are in agreement with the CHARA and AMBER 
measurements at long baselines, where all the circumstellar mat- 
ter appe^^Jo_bereTOlyed 2 and are similar to the values adopted 
bvrTannirkula m et alj d2008l) . The HD 163296 SED is shown in 
Figs.[8]and0 



4.2. incoherent flux 

In our models, we consider incoherent fl ux possibly em i tted b y 
an extended halo, as first suggested by Monn ieTet al.l (|2006). 
They estimated its contribution to be 5% of the total H band 
fluxes (i.e., 10% in V 2 ), a value that T08 also used when mod- 
elling their K band data. Our H band data suggest a slightly 
higher value of ~8% (i.e., ~15% in V 2 ). We adopt this value 
when fitting both the H and K band interferometric observations. 
The precise origin of this emission is unknown and discussing it 
is beyond the scope of this paper. 

4.3. Disk position angle and inclination 

Values of inclination (i) and position angle (PA) of the 
HD 163296 disk (i.e., its major axis) have been de- 
rived at different wavele ngths with a varie ty of tech- 
niques and ten d to agree: IWassell et ail d2006l) (i=51°+g'°; 
PA=139°+15°), llsellaet all d2007b (i=46°+4°;PA=128°±4°) 
and T08 (i=48°+2°;PA=136°+2°). They are also consistent with 
our visibility data, which infer an inclination of 40° + 10° and a 
position angle of 140° ±15°, when fitted with geometrical mod- 
els of uniform brightness ring. Since we do not constrain these 
parameters better, in the following, we adopt the values derived 
by T08, i.e., i=48° and PA=136°. In the following sections, we 
use the effective baseline, which is defined by : 



B e ff = B ^cos 2 (9) + cos 2 (®)sin z (6) 

where 9 is the angle between the baseline direction and the major 
axis of the disk, and <b is the disk inclination (T08). This repre- 
sentation allows us to show the data in a concise way once the 
inclination and position angle of the disk are known. 

4.4. A disk rim 

Early NIR interferometric studies of Herbig Ae stars have shown 
that standard accretion disks, extending u p to the dust subli- 
matio n radius do not fit the observations ( Millan-Gabet et all 
2007) and that superior fits are obtained by assuming that the 
disk develops a curved rim, probably (but not necessarily) con- 




Effective spatial frequency (m/|im) Effective spatial frequency (m/|j,m) 

Fig. 7. Model visibilities when considering a rim alone, located at 0.23, 0.36, and 0.45 AU (dashed red, dotted blue, and full black 
lines respectively), versus effective spatial frequency B e s/A, compared to the measured broad-band visibilities (blue full squares) 
in the K band (left) and the H band (right). The predictions of the model described in T08, which includes an additional inner 
component interpreted as gas, is given in red (T08; dash-dot-dot). The K band measurements obtained at CHARA and Keck-I are 
overplotted (circles and triangles respectively) in the left panel, as well as the IOTA H band data (triangles) in the right panel. An 
incoherent flux contributing to ~8% of the total H and K band fluxes is considered. 



trolle d by dust sublimation and its dependence o n gas den- 
sity dlsella & Nattal 120051: iTannirkulam et alj 120071) . Not only 
the disk rim hypothesis is supported by ph ysical calculations 
dDullemond et alj200lHfsella & Nattal2 005):' the overall proper- 
ties of HD 163296 are also consistent with this model. The SED 
appears to be consistent with a rim emis sion up to 7-8 nm and a 
disk in the shade at longer wavelengths dDullemond et al]|200lt 
llsella et alJ[2007h . In addition, our measured visibilities show a 
dependence on hour angle (or similarly on baseline position an- 
gle Bpa) consistent with an elliptical shape, as expected in either 
a rim or any axisymmetrical distribution seen at high inclina- 
tion. Finally, the detection of non-zero closure phases supports 
the idea that an asymmetric brightness distribution - such as the 
disk inner rim - contributes to the NIR emission. 
As in the NIR, th e contribution of the disk outside the rim 
can be neglected dlsella et alJ [2007b . we begin by examining 
a star + rim system to model the NIR emission following 
Isella & Natta (2005). We adopt a stellar luminosity and mass 
of 30 L and 2.3 M G , respectively, and an effective tempera- 
ture of 9250 K. The disk is assumed to be in hydrostatic equi- 
librium. The dust consists of silicate s with optical properties 
given by Weingartner & Draine (2001). The evaporation temper- 
ature, of around 150 K, depends on the local gas density as in 
Polla ck et al.l (fl994). Since the shape of the rim is controlled by 
the largest grains, we consider a single size for the silicate dust, 
which is therefore the only free parameter in the model. The de- 
pendence of the evaporation temperature on z implies that the 
distance from the star at which dust evaporates increases with z, 
i.e., that the rim is curved. 

As can clearly be seen in Fig. [7] our interferometric mea- 
surements are inconsistent with a circumstellar emission pro- 
duced by a rim only, regardless of its location. This figure shows 
the AMBER broad-band visibilities as a function of the effec- 
tive spatial frequency B e s/A. In addition to the AMBER data, 



we include the K band v isibilities obtained at C HARA (circles; 
I Tannirkulam et alj2 008) and at Keck-I (triangles; lMonnier et al.l 
2005 ) as well as the H band IOTA data (upside down triangles; 
Mon nier et al.f 2006). We emphasize that all measurements from 
these four interferometers are compatible with each other within 
their error bars. The results for three rim models computed by 
assuming a 30 L stellar luminosity and silicate grain sizes of 3, 
0.6, or 0.3 yt/m, corresponding to inner rim radii at 0.23, 0.36, and 
0.45 AU, respectively, are overplotted. These models produce an 
emission that does not exceed the observed flux at any wave- 
length, and, contributes respectively, to about 80, 50, and 28% of 
the observed emission in H, and 91, 86, and 66% in K. All mod- 
els show large departures from the observations. In addition, all 
these models have very asymmetric emission, and produce a clo- 
sure phase signal greater than observed (see Fig.Q~2l red dashed 
lines). Our large data sets, which include H and K visibilities 
and closure phases, reinforce the conclusions of other authors, 
namely that there should be an additional, symmetric emission 
closer to the star (T08). 

4.5. A bright inner disk 

A rim alone cannot reproduce our observations. As a matter of 
fact, the lack of bounces in the visibilities suggests a continu- 
ous and smooth distribution of matter in the inner regions of 
HD 163296, and a low contribution of the rim to the NIR emis- 
sion. In the following, we explore the effect of adding emission 
from an inner disk (inside the rim) to both the star and rim contri- 
butions. We fix the rim location to be Rri m =0.45 AU. This large 
radius is implied by the results of our three-component models, 
as we show in the following. It can, for example, be obtained by 
considering small silicate grains (from 0.05 to 0.5 fim) in a high- 
density dusty disk, or micron-sized grains in a low-density re- 
gion, where the evaporation temperature is lower. In both cases, 
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the rim effective temperature is quite low (T eff ~1000 K) and 
accounts for the total observed flux in the wavelength range 4- 
8 fim, but not at shorter wavelengths (see Fig. [8] cyan dashed 
line). The SED of the proposed additional inner disk can be com- 
puted by subtracting both the stellar and disk rim fluxes from the 
observed ones. As shown in Fig. [8] (dotted line), it peaks in the 
H band and is higher than the rim emission at all wavelengths 
shorter than about 2.5 fim. Therefore, the NIR emission is dom- 
inated by the inner disk component rather than the rim. 
We obtain a first estimate of the emission properties of the inner 
disk by modeling it as a region of constant surface brightness be- 
tween an inner and outer radius. The surface brightness is con- 
strained by the condition that the integrated flux must be equal to 
the flux derived from the SED for the additional component. We 
compute visibilities and find that these three-component models 
reproduce the H and K band observations quite well over the 
entire range of baselines if the smooth, inner emission extends 
between ~0.10 and ~0.45 AU. The emission can be roughly de- 
scribed as that of a diluted black-body with temperature ~ 1600 
K and optical depth ~0.2, which decreases as A~ l 6 with in- 
creasing wavelength. A natural interpretation is that the emis- 
sion originates in an optically thin region inside the rim. In the 
following sections, we discuss the physical nature of this com- 
ponent. 

4.5.1 . An inner gaseous disk ? 

Several authors have suggested that NIR interferometers de- 
tect the emission of gas within the dust evaporation radius 
dEisner et al] l2007h llsella et alj l2008t iTannirkulam et all l2008t 
iKraus et al.ll2008bl) . This interpretation presents several difficul- 
ties in the case of HD 163296, where the additional component 
dominates the emission in H and K, and where good observa- 
tions of high spectral resolution exist over a wide range of wave- 
lengths. 

Models of the emission of purely gaseou s disks inside 
the du st sublimation radius were computed by Muz erolle et alj 
(2004) for HAe stars assuming LTE opacities. For typical HAe 
accretion rates (~ 10~ 7 M yr _1 ), the inner, dust-free, disk gas 
surface density is higher than 10 3 g/cm 2 , the gas temperature 
ranges from a few thousand to a few hundred K, and the gas is 
fully molecular. The NIR emission is sufficiently strong to, in 
principle, account for the observations. However, as has been 
pointed out by several authors, the models also predict many 
strong molecular bands (mostly water and CO overtone transi- 
tions ) that are absent in the HP 163296 spectrum (see Fig. [8j 
left) (Naiita et al. 2000, 2007). This problem is also discussed in 
Naiit a et al.1 d2009l) . in connection with the non-detection of CO 
overtone and water emission in the Herbig Ae star, MWC480, 
that also exhibits a hot compact NIR excess previously de- 
tect ed with inte rferometry and interpreted as resulting from wa- 
ter (|Eisne3|2003). 

The assumption that the gas is in LTE, however, is certainly 
inappropriate at least in the upper disk layers, where the stel- 
lar radiation can penetrate, ionize, and dissociate matter well 
above the LTE predictions. However, these thin gas layers are 
unlikely to contribute significantly to the broadband observed 
fluxes. We computed the emission from thin layers of gas, us- 
ing the code Cloudy (Ferlan det al.ll 1998b . We assumed that the 
region inside the rim can be described as a geometrically thin 
disk, heated from the outside by the star. The disk extends from 
0. 1 to 0.45 AU and has a constant surface density. In these con- 
ditions, as long as the disk is optically thin to the stellar radia- 
tion, the gas is mostly atomic and the ionization fraction is low 



(e.g., < 0.01). In the NIR, H dominates the emission. We show 
in Fig. [8] right, that for surface density values of 0.1, 1, and 6 
g/cm 2 (all much lower than predicted by the viscous accreting 
disks modelled by Muzerolle et al. 2004), the continuum emis- 
sion is always too weak to be significant. Increasing it further 
would require higher surface densities, in which case LTE con- 
ditions would very likely be reached. In fact, for 2 = 6 g/cm 2 , 
the mean optical depth to the stellar radiation is already of or- 
der unity. Moreover, we note that the wavelength dependence of 
the non-LTE continuum is inconsistent with the observations. In 
particular, the H- bound-free emission produces a sharp drop at 
1 .6 fim corresponding to its activation energy , which would be 
seen in the HD 163296 high resolution spectra dSifko et al.l2 008) 
if the gas emission was higher. 

T08 suggested that hot gas is the physical origin of the addi- 
tional inner component. We note that their models cannot repro- 
duce our AMBER/VLTI data, mostly because of their small in- 
ner rim radius. However, we examined the possibility that the gas 
is heated to higher temperatures by additional energy sources. 
We computed the properties of the same disk model but after 
fixing the temperature to 8000 K. In this case, the ionization frac- 
tion is high (»0.5) and the emission is dominated by bound-free 
processes. The emission can be very high, but its wavelength de- 
pendence is inconsistent with the observations (see Fig.[8j right). 

Based on these crude considerations, we tend to exclude that 
the NIR flux detected by the interferometers is dominated by the 
emission of hot gas inside the dust sublimation radius. However, 
it is clear that, before this can be definitely ruled out, one needs 
more realistic, non-LTE models that treat the transition from op- 
tically thin to optically thick layers, i.e., from atomic to molec- 
ular gas in a dust-free environment. These models would also 
be important for the interpretation of the hydrogen and helium 
recombination lines, which appear very strong in the models. 
Atomic lines - mostly hydrogen and helium ones - are often in- 
terpreted as being emitted in magnetospheric accretion columns 
of gas. This, however, is unlikely to be true for most Herbig Ae 
stars, based on the results obta ined with spectro-interferometry 
aroun d the Br r emission line dKraus et alj|2008al : lEisner et al.l 
2009), since this line seems to be formed in most cases, in disk 
material closer to the star than the silicate dust sublimation ra- 
dius but outside the corotation radius. 

4.5.2. An inner dusty disk ? 

Inside the silicate sublimation radius, not only gas but also more 
refractory grains can exist and contribute to the observed NIR 
emission. In this section, we explore the possibility that a layer 
of refractory grains, extending between an inner and an outer 
radius (inside the rim), accounts for the interferometric and pho- 
tometric observations of HD 163296. We assume that the layer 
is optically thin in the vertical direction with surface density pro- 
file Sdust and that its vertical optical depth is proportional to l/r, 
where r is the distance from the star (D'Ales sio et al . 1999). We 
compute grain temperatures and emissivity in the H and K bands, 
and vary the inner and outer radii (R; n and R out , respectively) 
as well as the optical depth of the layer. We then computed the 
emission, visibilities, and closure phases for models that include 
this inner layer, the star, and the rim. We consider separately 
three grain species known to be refractory: iron, graphite, and 
corundum (aluminium oxides). The iron and corundum grain 
cross-sections are computed for spherical grains f rom the op- 
tical c onsta nts tabulated in iPollack et alj dl994l) . iKoike et all 
d 19951) . and Begem ann et alj dl997l) . We use the graphite cross 
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Fig. 8. SED of HD 163296 dSitkoet all 120081: iTannirkulam et alj|2008h with the predictions of gaseous disk models. The flux 
emitted by the rim located at 0.45 AU (cyan dashed line) as well as the additional inner disk emission (black dotted line) are added 
in bot h panels to allow direct comparisons. Left: the predictions for optically thick and dense layers of gas in LTE dMuzerolle et alj 
2004) are overplotted (red full line). In this case, the gas is mostly in a molecular state. Right: the continuum emission as predicted 
by thin disks of gas in non-LTE, heated from the top by the stellar radiation, are shown. The gas extends from 0. 1 to 0.45 AU and 
has constant surface densities of 0.1, 1, and 6 g/cm 2 (red line with squares, blue line with triangles and green line with diamonds, 
respectively). In these conditions, the gas is mostly atomic. The predictions of a fully ionized layer of gas at 8000 K with a constant 
surface density of 0.06 g/cm 2 are also reported (magenta curve with diamonds). As long as the disk remains optically thin, the 
emission is proportional to the surface density, and this value was chosen only to display its wavelength dependence. 



se ctions tabulated by B. Drainfl based on the optical constants 
of lLaor & Draind (Il993l) . For relatively small grains, the radial 
temperature profile follows oc r~° 4 (Fig. [9] left). For comparison, 
we also show the temperature profile of very large grains (which 
behave as blackbodies), which have flat opacity from the UV to 
the IR (T oc r -0 - 5 ). 

These models are very simple, but probably not unrealistic. The 
strongest approximation concerns the dust temperature, which 
we compute by assuming that each grain is heated by the stel- 
lar radiation, and attenuated by an average optical depth T mr oc 
0.25. In fact, the attenuation is not constant, since the optical 
path toward any individual grain varies not only with r but also 
with the incident angle of the stellar radiation. Once the tem- 
perature of the grains is known, then the emission is computed 
at all NIR wavelengths using a ray-tracing algorithm. As dis- 
cussed in the following, a rim caused by silicate condensation 
can form in the low density region that we propose and its prop- 
erties may also affected by refractory grains in the inner disk, 
which absorb a small fraction of the stellar radiation. However, 
a self-consistent calculation of the rim properties is beyond the 
purpose of this paper. In t his section, we mod el the rim following 
llsella & Nattal d2005l) and llsella et alj (120061) . assuming micron- 
size silicates, an evaporation temperature of -1250 K, and an 
effective stellar luminosity of 75% L+ to account for the shield- 
ing by the inner disk. The rim radius is about 0.45 AU and its 
effective temperature is about 1000 K. 

We obtain a good fit to the SED, as shown in Fig. [9] (right), 
for R in = 0.10 AU and R om = 0.45 AU, and an optical depth in 
the H band at of 0.31, 0.2, and 0.25 for iron, graphite, and 



2 http://www.astro.princeton.edu/ draine/dust/dust.diel.html 



corundum, respectively. Table |4] indicates the grains properties 
used in these models. The total mass within 0.45 AU, which is 
always very low (< 10~ 6 M ffi ), is also reported. Although we 
have not tried to constrain the grain parameters in any detail, we 
note that very large refractory grains do not provide an equally 
good fit to the SED, as the ratio of the H to the K emission is 
always too high. 

These models successfully reproduce the interferometric data, in 
terms of both visibilities and CR Since all three models produce 
similar results, we present the case for iron grains. Figure [10] 
shows the broad-band visibilities compared with the predic- 
tions of such a model. Figure Q~2] presents the broadband clo- 
sure phases plotted versus hour angle and the model predictions 
calculated for each telescope configuration. We overplot the pre- 
dictions of the rim-alone models (red, dashed line) to show how 
adding the optically thin dusty, inner disk emission smoothes out 
the asymmetry induced by the rim by the correct amount. Table|5] 
summarizes the structural parameters of this model, as inferred 
from SED and interferometric data fitting, and Fig.QT]shows the 
corresponding H and K band images. 

The gas density of the inner disk can be derived from the 
values of Sdust, once the abundances of iron, carbon, and alu- 
minium in the solid species are known. Assuming, for example, 
that 50% of the iron is in grains, the gas surface density at 
0.10 AU will be 0.2 g/cm 2 (0.02 g/cm 2 if 30% of carbon is in 
graphite, or 1.3 g/cm 2 if all aluminium is in corundum). The 
gas density can be higher if a lower fraction of the metals is 
condensed. However, it seems likely, from the considerations of 
Sect. 4.5.1, that the gas density cannot be too high. A density 
only a few times higher than the above lower limits is inferred 
by the properties of the rim in HD 163296. In particular, both 
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Fig. 9. Left: radial temperature profiles for graphite (green line with triangles), corundum (red line with crosses), or iron (blue line 
with squares). The prediction from a black body (i.e., large grains) is overplotted to allow comparisons. The vertical dotted line 
indicates 0. 10 AU. Right: SED of HD 163296 together with the predictions of three models that include an inner disk made of dust 
in addition to the star and the rim. The additional disk component has structural characteristics that are reported in Table [5] and 
is made of a single refractory species - either graphite, corundum, or iron. The dotted lines represent the corresponding emission, 
while the full lines show the total flux. 




1.0 



> 



0.5 



0.0 





i i i i i 

H band 




v IOTA 




















I.I.I. — 1 



20 



40 



60 



Effective spatial frequency (m/|j,m) Effective spatial frequency (m/|j.m) 

Fig. 10. Left: Visibilities produced by our model that assumes that iron grains are responsible for the additional dusty inner disk 
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the large rim radius and its low effective temperature can be 
reproduced in a low density dis k by silicates of micron size, 
as inferred in several HAe stars (llsella et al.ll2006l) . and do not 
require very small grains. Assuming that the rim is produced by 
the evaporation of silicates and that all the silicon is in olivine 



of ~ lyum size, we analyzed the rim properties in disks of 
increasing (but still low) surface density. For 2 = 5(r/0.1AU)~ l 
g/cm 2 , the vertical optical depth to the stellar radiation is ~10 
for a gas surface density of 1 g/cm 2 at R T i m , large enough 
to allow the formation of an optically thick rim that can be 
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Fig. 11. Model images of the circumstellar matter surrounding HD 163296, in the case where the inner disk is consists of refractory 
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Table 4. Refractory dust model parameters. For each species, 
the optical depths in H and K bands at 0. 10 AU (R; n ) are reported 
as well as the minimum and maximum grain sizes considered in 
the calculations. The resulting surface density £dust and the total 
mass of dust grains are also given. 
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[pm] 
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8.7 10~ 8 


Corundum 
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0.25 
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3.2 10~ 4 


1.3 10~ 6 



Table 5. Parameters of our model with iron grains distributed 
between R; n and R r ; m , the location of the disk rim. For both 
bands, the ratio of the stellar, rim, and inner disk contributions 
(F„, F r ; m , and F respectively) to the total flux in the model (F tot ) 
are reported. 
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modelled following llsella & Nattal (2005). The low gas density, 
and the correspondingly low evaporation temperature, moves 
the rim radius further from the star. Assuming a reasonable scale 
height of 1CT 3 - 1CT 2 AU, the gas density is ~ 10~ n g/cm 3 , and 
the silicate evaporation temperature of the order of 1 150 K. The 
corresponding rim radius is ~0.4-0.5 AU, and its effective tem- 
perature is about 1000-1 100 K, as required to fit the HD 163296 
SED and interferometric data. 

The only difficulty in assuming that the inner disk emis- 
sion originates in grains within the silicate sublimation radius, 
is the need for them to survive at temperatures much higher 
than is generally assumed. The three types of grains that we 



examined reach temperatures of 2100-2300 K at ~ 0.10 AU 
(Fig. [9l left). While similar values are possible for graphite 
(Krugel 2003), they are too hig h for both iron and corundum in 
the p r essure of our inne r disk dPollack et alJll994l; iPosch et all 
2003; Kama et al. || 2009l) . However, there is room for discussion 
dNajita et alj|2009t) . as the balance between gas and dust in the 
conditions of the i nner disk should be reconsidered in detail 
dDuschl et alj[l9 99). We emphasize that only a small amount of 
refractory grains need to survive these high temperatures, i.e., 
probably a minor fraction of the original population. The grains 
that we have considered are likely candidates, but if other, more 
refractory species can form, they would certainly fit the observa- 
tions equally well. 



5. Summary and conclusions 

This paper had discussed the largest set of NIR interferometric 
data collected so far for a young star. HD 163296 has a well 
studied disk at large spatial scale, which motivated our interpre- 
tation of the NIR interferometry using a star + inner disk model. 
Both interferometric and photometric data can be accounted for 
by an inclined disk with a low density inner region. The NIR 
continuum emission is then not dominated by the thermal emis- 
sion from the dusty disk rim located at the sublimation radius of 
astronomical silicates, but by the additional optically thin com- 
ponent located inside. This component emits about 32% of the 
stellar luminosity, but 54% (50%) of the observed radiation in 
H (K) band. Because of the silicate condensation and the strong 
increase in opacity, the disk rim forms in this low-density region 
at ~ 0.45 AU from the star and emits about 16% (36%) of the 
H (K) band flux. The combination of the unresolved stellar ra- 
diation, the smooth and point-symmetric emission of the inner 
disk region, and the skewed disk rim emission can successfully 
explain the visibilities and non-zero closure phases measured in 
both bands. 

The nature of the emission in the inner disk remains a mat- 
ter of discussion. We argue against gas being mainly responsible 
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Fig. 12. Closure phases versus hour angle. In each panel, the predictions of the disk rim model (dashed red lines) and our model that 
includes iron grains (full black line) are overplotted. The corresponding telescope configuration and wavelength band are indicated 
in the upper, left corners. The K band and H band measurements are plotted with circles and triangles, respectively. 



for this continuum emission. A dense, cold disk, as expected for 
viscous accretion models (Muzerolle et al. 2004), would pro- 
duce strong molecular lines that are not seen in high-resolution 
spectra. Non-LTE tenuous gas layers, in an atomic state if only 
heating by the star is assumed, or fully ionized by additional 
sources of energy, cannot account for the observed properties of 
the NIR continuum. However, we emphasize that self-consistent 
models of dust-free gaseous disks are not currently available, but 
are needed to exploit the full potential of the interferometric ob- 
servations. 



We suggest instead that a small fraction of refractory grains 
survive very close to the star. We propose models for the opti- 
cally thin emission of the innermost region, using various kinds 
of refractory grains, distributed from 0.10 AU to 0.45 AU. The 
dust surface density provides only a lower limit to the gas sur- 
face density, as we do know neither the exact nature of the grains 
nor their abundance. However, we find that a low density region 
is consistent with the location and properties of the rim, as con- 
densation of silicates will occur, and with the lack of molecu- 
lar features in the spectrum of HD 163296. We expect the gas 
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in the inner disk indeed to be mostly atomic, in non-LTE, and and S. Renard for fruitful discussions. We acknowledge the anonymous referee 

although its continuum emission will be weak, hydrogen lines for his comments that improved the clarity of the paper, 
can be strong. The models used in Sect. 4.5.2 to argue for the 
presence of refractory grains, make a number of crude assump- 
tions, and improved models that self-consistently compute the 
grain temperature and emission in the thin disk as well as the 
properties of the rim are being developed. Our study indicates 
that the inner region of HD 163296 is quite empty, with a very 
low surface density that is inconsistent with a dense accretion 
disk. For comparison, a surface density of 1 g/cm 2 at ~0.10 AU 
corresponds, in a standard accretion disk (a=0.01), to an accre- 
tion rate of ~ 10" 11 M^yr" 1 , much lower tha n typical values 
for Herbig Ae stars (Garcia Lope z et all 12006). With our data, 
we do not constrain the outer radius of this low density region, 
which can be larger than 0.5 AU. However, we know that at large 
radii the HD 163296 disk is massive and de nse, as shown b y the 



seems likely that HD 163296 has a dense disk with an inner cav- 
ity, and that we observe it just before it re aches the transition 
disk phase, as suggested bv lSitkoetal] (120081) . 

It is fair to emphasize that we make no claim that our 
interpretation is unique. As stated at the beginning of Sect. 4, 
we assumed that the NIR emission of HD 163296, at spatial 
scales of less than 0.5 AU, is dominated by the emission of a 
circumstellar disk. Moreover, we interpreted the non-zero values 
of the CP as evidence of the asymmetric emission of a disk 
rim. While the properties that we derived for the smooth inner 
emission are probably robust, the existence of the rim is less 
so. In particular, the lack of visibility bounces at large baselines 
argues against the presence of a rim. In this case, the observed 
closure phases may possibly be caused by any asymmetric 
brightness distribution, such as a symmetric flared disk with a 
stellar contribution that is off-centered by a few percent of th e 
inner disk radius with respect to the disk ( Ma lbet et al.l 1200 1 ). 
by a hot spot on the disk, or by a density discontinuity. Only a 
larger (u,v) coverage, providing access to more details of the 
morphology could solve this ambiguity. 

Our interpretation of the smooth, inner emission as originat- 
ing in refractory grains requires their survival at very high tem- 
peratures (~ 2100 - 2300 K), much higher than expected, even 
for the most refractory grains, at the pressure of the low-density 
inner disk dPollack et al.l Il994t iKama et all 120091) . However, 
clearing the innermost regions and optically thin emission from 
left-over refractory grains could be a common phenomenon in 
Herbig Ae stars. Their presence within the first few tenths of an 
AU is a promising interpretation of the observed depletions in re- 
fractory dust species - such as iron - in jets of young stars that ar e 
launched from this region dNisini et al.ll2"00H " lPodio et"aIll2006T) . 
Similar interferometric studies, with a large number of measure- 
ments in various wavelength bands simultaneously, should be 
performed for a large sample of stars. A higher level of com- 
plexity in models is also needed to account for both the dust and 
the gas emission in a self-consistent way. Finally, the advent of 
the next-generation of imaging instruments will hopefully pro- 
vide unambiguous constraints on these complex environments. 
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